m = 20;   % Image size m x m.
Rd = 3;   % Disk radius.
Center = (m+1)/2;
u0 = zeros(m,m);
for i = 1:m
    for j = 1:m
        if (Center - i)^2 + (Center - j)^2 < ExParam^2
            u0(i,j) = 1;
        end;
    end;
end;
for diskNum = 1:N
    DiskX(diskNum) = (m-1)*rand+1;
    DiskY(diskNum) = (m-1)*rand+1;
    u = u0;
    for i = 1:m
        for j = 1:m
            if (DiskY(diskNum) - i)^2 + (DiskX(diskNum) - j)^2 < Rd^2
                u(i,j) = 1;
            end;
        end;
    end;
    X(diskNum,:) = reshape(u,1,m*m);
end;
% Since this is a special manifold, plot separately.
axes(handles.maniAXES);
t = 0:0.1:2*pi+0.1;
plot(ExParam*cos(t)+Center,ExParam*sin(t)+Center);
axis([0.5 m+0.5 0.5 m+0.5]);
hold on;
handles.ColorVector = (sqrt((DiskX-Center).^2+(DiskY-Center).^2))';
scatter(DiskX,DiskY,12,handles.ColorVector');
hold off;
handles.X = X;
updateString{1} = 'Occluded Disks example loaded.';
updateString{3} = 'Warning: Embedding may be slow.';
